On the intermittency exponent of the turbulent energy cascade 
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We consider the turbulent energy dissipation from one-dimensional records in experiments using 
air and gaseous helium at cryogenic temperatures, and obtain the intermittency exponent via the 
two-point correlation function of the energy dissipation. The air data are obtained in a number of 
flows in a wind tunnel and the atmospheric boundary layer at a height of about 35 m above the 
ground. The helium data correspond to the centerline of a jet exhausting into a container. The air 
data on the intermittency exponent are consistent with each other and with a trend that increases 
with the Taylor microscale Reynolds number, R\ , of up to about 1000 and saturates thereafter. On 
the other hand, the helium data cluster around a constant value at nearly all R\, this being about 
half of the asymptotic value for the air data. Some possible explanation is offered for this anomaly. 

PACS numbers: 47.27.Eq, 47.27.Jv, 47.53.+n, 05.40.-a, 02.50.Sk 



I. INTRODUCTION 

That turbulent energy dissipation is intermittent in 
space has been known since the seminal work of Batchc- 
lor & Townsend [jj . The characteristics of intermittency 
are best expressed, at present, in terms of multifractals 
and the multiplicity of scaling exponents; see e.g. Ref. Q. 
In the hierarchy of the scaling exponents, the so-called 
intermittency exponent characterizing the second-order 
behavior of the energy dissipation, is the most basic. In 
this paper, we address the Reynolds number variation of 
the intermittency exponent and its asymptotic value (if 
one exists). 

The intermittency exponent has been determined by 
a number of authors in the past using several different 
methods; for a summary as of some years ago, see Refs. 
0,0. Recently, Cleve, Greiner & Sreenivasan Q evalu- 
ated these methods and showed that the best procedure is 
to examine the scaling of the two-point correlation func- 
tion (e(x + d)e(x)) of the energy dissipation e. Other 
procedures based on moments (ef) and (ln 2 £;) — (lne;) 2 
of the coarse-grained dissipation Ei — (1/7) J l e(x)dx, or 
the power-spectrum P(k) — \e(k)\ 2 , are corrupted by the 
unavoidable surrogacy of the observed energy dissipation. 

The follow-up effort [|| was able to characterize and un- 
derstand the functional form of the two-point correlation 
function beyond the power-law scaling range. Within 
the theory of (binary) random multiplicative cascade pro- 
cesses, the finite-size parametrization 



ilar to the integral scale (as discussed later) but typi- 
cally larger. Here, \x is the intermittency exponent, and 
the first term is the pure power-law part and the sec- 
ond term is the finite-size correction. Comments on the 
constant c will be given further below. The one atmo- 
spheric and the three wind tunnel records employed in 
Ref. @ were found to be in accordance with this close- 
to-universal finite-size parametrization. Even for flows at 
moderate Taylor-microscale Reynolds numbers R\, Eq. 

allowed an unambiguous extraction of \x. For the 
four data records, a weak dependence of /i on R\ was 
observed in but was not commented upon in any de- 
tail. The present analysis examines many more records 
and attempts to put that Reynolds number dependence 
on firmer footing. 



(e(x + d)e(x)) 



(1) 



was derived, introducing the cascade length L casc as a 
meaningful physical upper length scale, this being sim- 



II. THE DATA 



Two of the data records examined here come from the 
atmospheric boundary layer (records al and a2) |7j, eight 
records from a wind tunnel shear flow (records wl to w8) 
@ and eleven records from a gaseous helium jet flow 
(records hi to hll) 0. We find that all air experiments 
show an increasing trend of the exponent towards 0.2 
as the Reynolds number increases. The helium data, 
on the other hand, show the exceptional behavior of a 
Reynolds- number- independent "constant" of around 0.1. 
Some comments on this behavior are made. 

It is useful to summarize the standard analysis in order 
to emphasize the quality of the data. As is the standard 
practice, energy dissipation is constructed via its one- 
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dimensional surrogate 

e(*) = 15,(J)\ (2) 

where the coordinate system is defined such that u is the 
component of the turbulent velocity pointing in the longi- 
tudinal x direction (the direction of mean motion). The 
coefficient v is the kinematic viscosity. Char acter istic 
quantities such as the Reynolds number R\ = ^/ {u 2 )X/v, 
based on the Taylor microscale A = \J (u 2 ) / {{du/dx) 2 ) , 
the integral length L, the record length L r ecord) the res- 
olution length Ax and the hot-wire length l w in units 
of the Kolmogorov dissipation scale 77 = (^ 3 /(e)) 1 / 4 are 
listed in Table [I] To calculate the numerical value of A 
the method described in 0] has been used. The inte- 
gral length is defined as the integral over the velocity 
autocorrelation function (using Taylor's hypothesis). In 
the atmosphere, where the data do not converge for very 
large values of the time lag, the autocorrelation function 
is smoothly extrapolated to zero and the integral is evalu- 
ated. This smoothing operation towards the tail does not 
introduce measurable uncertainty in L. The energy spec- 
trum, which is illustrated in Fig. ^ for the records a2, wl 
and h7 as representatives of the three different flow ge- 
ometries with Reynolds numbers ranging from the small 
to the very large side, follows an approximate — 5/3-rds 
power over the inertial range. The wind tunnel records 
are relatively noise-free, while the helium data are af- 
fected by instrumental noise significantly, as evidenced by 
the flattening of the energy spectrum for high wave num- 
bers. The atmospheric data fall somewhere in-between. 
The effect of this high-frequency noise, and of removing 
it by a suitable filtering scheme, will be discussed later. 
Again for the three representative cases a2, wl and h7, 
Fig. [3 shows the two-point correlation defined in Q for 
the surrogate quantity J5J and compares it to the best-fits 
for the proposed form of finite-size parametrization given 
by Eq. As is also the case for all the other records, 
the agreement is quite substantial and unambiguous. The 
upturn at small separation distances d < A* has been 
explained in [j| as the effect of the surrogacy of the en- 
ergy dissipation. Hence, we call A* the surrogacy cutoff 
length. The two-point function decorrelates at a length 
scale L casc that is substantially larger than the integral 
length scale. The cascade length L casc and the surro- 
gacy cutoff length A* are also listed in Table [I] for all the 
records inspected. 



III. INTERMITTENCY EXPONENT 

To determine the intermittency exponent /x, the data 
are fitted to Eq. Q using a best-fit algorithm (see Fig. 01 . 
These values are also listed in Table H] The local slopes 
from the best fits, plotted as insets in the figure, show 
that deviations from the pure power-law become evident 
only towards large values of the separation distance d. 




f(s" 1 ) 



FIG. 1: Power spectral densities of the velocity fluctuations 
for the records a2, wl and h7. The sharp peaks are artifacts. 
For the record h7, the filtered power spectral density is also 
shown. 
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TABLE I: Taylor-microscale based Reynolds number Rx, the integral length scale L in units of the Kolmogorov scale 77, the 
record length I/ rccor d, the Taylor microscale A, the resolution scale Ax (= sampling time interval x mean velocity), the length 
of the hot wire l w , the intermittency exponent /x, the cascade length Lease and the surrogacy cutoff length A* for the two 
atmospheric boundary layer data (al,a2) [3, the eight wind tunnel data (wl-w8) y| and the eleven sets from gaseous helium 
jet (hl-hll) measurements. 



data set 


Rx 


L/r, 


-^record / L 


A/7? 


Ax/7/ 


lw/v 


H 


Lease 1 f] 


A*A? 


al 


9000 


5xl0 4 


1000 


187 


1.29 


1.755 


0.216 


322743 


3.9 


a2 


17000 


7.5 xlO 4 


970 


246 


3.64 


1.534 


0.202 


509354 


9.1 


wl 


208 


539 


28000 


27 


2.42 


1.052 


0.143 


1164 


24 


w2 


306 


484 


102500 


35 


1.98 


1.780 


0.155 


1939 


26 


w3 


410 


697 


127700 


38 


2.67 


2.533 


0.151 


2707 


27 


w4 


493 


968 


193500 


44 


2.79 


3.382 


0.145 


3228 


31 


w5 


584 


1095 


88600 


44 


2.71 


0.890 


0.172 


4062 


27 


w6 


704 


1365 


117700 


48 


2.90 


1.079 


0.173 


4343 


29 


w7 


860 


1959 


89500 


53 


2.63 


1.580 


0.176 


5513 


26 


w8 


1045 


2564 


77500 


64 


2.97 


1.927 


0.171 


7469 


27 


hi 


85 


102 


197000 


22 


1.20 


0.040 


0.12 


934 


10.8 


h2 


89 


101 


175000 


22 


1.05 


0.025 


0.128 


472 


10.5 


h3 


124 


165 


100000 


26 


0.98 


0.068 


0.102 


738 


17.7 


h4 


208 


344 


85200 


33 


1.75 


0.088 


0.154 


1258 


14.0 


h5 


209 


277 


59000 


23 


0.97 


0.072 


0.083 


1559 


10.7 


h6 


352 


606 


62400 


47 


2.25 


0.165 


0.13 


2254 


22.5 


h7 


463 


1011 


32100 


50 


1.93 


0.310 


0.092 


10438 


25.1 


h8 


885 


1442 


40100 


47 


3.45 


0.763 


0.089 


18954 


10.3 


h9 


929 


2064 


29800 


48 


3.67 


0.762 


0.079 


8434 


18.3 


hlO 


985 


2144 


37800 


48 


4.83 


0.837 


0.105 


23659 


14.5 


hll 


1181 


3106 


26900 


57 


4.97 


1.097 


0.061 


14921 


19.9 
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FIG. 2: Best fits of expression Q to two-point correlators 
extracted from data sets a2, wl, unfiltered h7 and filtered h7 
(h7f). The inset figures show the local slopes, compared with 
the fits given by (1). For comparison, power-law fits with 
the extracted intermittency exponent listed in Table I (and 
fi = 0.13 for h7f) are shown as dashed straight lines drawn 
with arbitrary shifts. 



The dashed lines in each figure are pure power-laws (with 
arbitrary shift) for the values of /j listed in the Table. 

For the atmospheric boundary layer, the analysis of 
two data sets yields a value of about 0.2 for the intermit- 
tency exponent /i; see Fig.|3J Note that, since the data set 
al contains both the longitudinal and transverse velocity 
components, one can form different forms of the surro- 
gate energy dissipation. It was found in Ref. 5| that all 
of them lead to the same value of the intermittency expo- 
nent. The filtering of the data has no measurable effect 
on the numerical value of /x. 

The wind-tunnel data wl-w8 seem to suggest a 
Reynolds number dependence of the intermittency ex- 
ponent for R\ of up to about 1000. The value /i = 0.2 
of the atmospheric boundary layer is reached only for 
higher R\. Unfortunately there is no laboratory data 
for R\ > 1000 so that there is a gap between the wind 
tunnel data and the atmospheric boundary layer data. 
Nevertheless, all the air data taken together appear to 
be consistent with a trend that increases with the Taylor 
microscale Reynolds number up to an R\ of about 1000, 
and saturates thereafter. This trend is also supported by 
results quoted in the literature H H GH HI 111 111 EM, 
although the finite-size form has not been employed 
for the extraction of the intermittency exponent. The 
literature values, shown in Fig. PJ fill the gap between 
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FIG. 3: The intermittency exponent, /i, extracted from a best 
fit of expression Q to the two-point correlator of the various 
data records, as a function of the Taylor-microscale based 
Reynolds number. Also shown are some values quoted in the 
literature. For some of the helium data, the lines show the 
shift resulting from the application of the Wiener filter to 
remove high-frequency noise. 



the present wind tunnel and atmospheric data. 

In contrast to the air data, the gaseous helium records 
hl-hll show a different behavior (Fig. It appears 
that, unlike the air data which show a gradual trend with 
R\, leading to a saturation for R\ > 1000, the helium 
data yield an intermittency exponent that is flat with 
R\ at a lower value of 0.1. It remains an open question 
as to why this is so. It would be important to settle 
this puzzle and clarify if this special behavior has other 
consequences for the helium jet data. 

To make some progress, we examined the helium data 
more closely. Perhaps the instrumental noise, which is 
seen in Fig. by the flattening of the energy spectrum 
for high frequencies, affects the accuracy of the calcula- 
tion of the energy dissipation. To account for such ef- 
fects, we applied a Wiener filter to the data, see again 
Fig. and recomputed the two-point correlation; the 
result is shown in Fig. [2Ji- The quality of the agreement 
with the finite-size parametrization remained the same 
but the numerical value for the intermittency exponent 
altered. Filtering produces different amounts of shift for 
different sets of data; see again Fig. [3] The most extreme 
change of the numerical value was found for h7, where the 
intermittency exponent changed from /i = 0.09 in the un- 
filtered case to \x = 0.13 in the filtered case. The differ- 
ence between the two values can perhaps be taken as the 
bounds for the error in the determination of the inter- 
mittency exponent. Given this uncertainty, one cannot 
attribute any trend with respect to the Reynolds num- 
ber for the helium data, and an average constant value of 
fi w 0.1 seems to be a good estimate for all helium data. 

Further questions relate to the spatial and temporal 



resolutions of the hot wire. The temporal resolution in 
the helium case is comparable to that in the air data (see 
Table [J; and, if anything, the ratio of the wire length to 
the smallest flow scale, namely l w /"q, is better for helium 
experiments. However, an important difference between 
the air data and the helium data concerns the length to 
the diameter of the hot wire. For air measurements, the 
ratio is usually of the order of a hundred (about 140 for 
al and a2 and about 200 for wl to w7), while it is about 
1/3 for hl-hll. In general, this is some cause for concern 
because the conduction losses from the sensitive element 
to the sides will be governed partly by this ratio, but the 
precise effect depends on the conductivity of the mate- 
rial with which the hot wires are made. For hot wires 
used in air measurements, the material is a platinum- 
rhodium alloy, while for those used in helium, the wire is 
made of Au-Ge sputtered on a fiber glass. This issue has 
been discussed at some length for similarly constructed 
hot wires of Emsellem et al. |T||. The conclusion there 
has not been definitive, but the helium data discussed in 
show another unusual behavior: unlike the air data 
collected in 0], the flatness of the velocity derivative 
shows a non-monotonic behavior with R\. See also fig- 
ure 4 of Ref. 0|. Whether the two unusual behaviors 
of the helium data are related, and whether they are in 
fact due to end losses, remains unclear and cannot be 
confirmed without further study. A further comment is 
offered towards the end of the paper. 

The data on the surrogacy cutoff length A* does not 
show a clear Reynolds number dependence. Referring 
again to Table [I] it appears that A* is directly related to 
neither A nor r\. However, definitive statements cannot 
be made because of the practical difficulty of locating A* 
precisely. 

Figure 0| illustrates the findings on the cascade length 
ratio L casc /r]. The ratio increases with the large-scale 
Reynolds number R = u' L caac /v as a power-law with 
the exponent of 3/4, exactly as anticipated if L C3SC were 
proportional to the integral scale. The ratio L casc /L is 
not exactly a constant for all the data (as can be seen 
from Tabic 1), but given the uncertainty in determining 
L and the absence of any systematic trend suggests that 
our supposition is reasonable. This is further reinforced 
by the variation of L casc /r] with respect to R\ (see inset), 
which also follows the expected behavior. We should note 
that it is difficult to single out the helium data in this 
respect. 

There is more to learn from fitting the finite-size 
parametrization Q to the experimental two-point cor- 
relations than merely extracting the intermittency expo- 
nent and the cascade length. As revealed by a closer 
inspection of the asymptotic behavior as d — > L casc , see 
again Fig. [2 the two-point correlation of the atmospheric 
boundary layer and wind tunnel records approach their 
asymptotic value of unity from above, whereas for most 
of the gaseous helium jet records the curve first swings 
a little below unity before approaching the asymptotic 
value. The expression is flexible enough to reproduce 



5 




flux density fed into the inertial-range cascade at the up- 
per length scale is somewhat larger than for the helium 
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FIG. 4: The dependence on the Reynolds number of the ra- 
tio of the cascade length Lease to the Kolmogorov scale r\. 
The a-, w- and h-records are represented by diamonds, circles 
and triangles, respectively. In the main graph, the straight 
line indicates a power-law scaling with exponent 3/4, and the 
Reynolds number is defined as R = u'L CB , sc /p. The inset 
shows the same data over R\ with the straight line indicating 
the expected power-law scaling of 3/2. The prefactor for the 
main graph is A — 1.3, for the inset A = 0.41. 



even this behavior. The derivation of within the the- 
ory of binary random multiplicative cascade processes, 
which has been presented in Ref. [g, also specifies the 
parameter 



(Tl 2 )(qLq R ) 
2<«£/n> 



- 1 



(3) 



in terms of cascade quantities. Normalized to (Id) = 1, II 
represents the initial energy flux density, which is fed into 
the cascade at the initial length scale L casc . (g^/ R ) and 
(qlQr) are second-order moments of the bivariate prob- 
abilistic cascade generator p(qL,lR) = p(qr,<1l), which 
we assume to be symmetric. Again the normalization of 
the left and right random multiplicative weights is such 
that (ql/r) = 1< Note, that log 2 (qf j , R ) = fi is equal 
to the intermittency exponent. Figure [5] shows various 
graphs of the two-point correlation with the expres- 
sion ||3J|, where parameters \i and L casc have been kept 
fixed, but c has been varied in the range < c < 1. 
We observe that for large c the two-point correlation ap- 
proaches its asymptotic value from above, whereas for 
small c it swings below one before it reaches the asymp- 
totic value from below. The transition between these 
two behaviors occurs at c w ti— ■ This translates to 
(U 2 )(q L q R ) = (2 1+ ^-l)/(l+/i), which is 1.08 for fi = 0.2 
and 1.04 for /i = 0.1. Hence, we are tempted to conclude 
that for the air data the fluctuation of the initial energy 




FIG. 5: Two-point correlator Q for various parameter values 
c = 0.6, 0.7, 0.8, 1/(1 + fi), 0.9, 1 (from left to right). The other 
parameters are fixed to Lease A? = 10 4 and /i = 0.20. 



jet data; this appears to be plausible and is one differ- 
ence between air and helium data. We also read this as 
an indication that {qbqR.) < 1, which is fulfilled if the 
left and right multiplicative weight are anticorrelated to 
some extent. As has already been discussed in a different 
context this anticorrelation is a clear signature that 
the three-dimensional turbulent energy cascade conserves 
energy. 



IV. CONCLUDING REMARKS 

In summary, we state that the picture of the turbu- 
lent energy cascade is robust and again confirmed by 
the excellent agreement between the two-point corre- 
lation density predicted by random multiplicative cas- 
cade theory and that extracted from various experimen- 
tal records. The cascade mechanism appears to be uni- 
versal, although its strength, as represented by the in- 
termittency exponent, seems to depend on the Reynolds 
number except when it is very high. The discrepancy 
between the air data on the one hand and the gaseous 
helium data on the other remains a puzzle (despite some 
possible explanations offered), and is in need of a fuller 
explanation. 
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